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ABSTRACT 


In  Technical  Report  4,^^  several  definitions  and  theorems  were 
presented  to  aid  in  the  analysis  of  a  sub-class  of  non-stationary  pro 
cesses  consisting  of  a  random  component  and  a  continuous  non-random 
function  of  time,  each  component  defined  over  the  same  finite  time 
interval.  This  report  supplements  the  results  of  Technical  Report  4 
both  by  extending  the  sub-class  of  non-stationary  processes  under  con 
sideration  and  by  including  additional  theorems. 
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References  ere  J'Ated  Mt  the  end  of  this  report. 
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I  INTRODUCTION 


This  report  is  divided  into  three  sections.  Section  I  summarizes 
some  of  the  results  of  Technical  Report  4.^  Section  II  describes  both 
theorems  that  are  natural  extensions  to  those  of  Technical  Report  4,  and 
extensions  to  the  sub-class  of  processes  considered  in  the  above  report. 
Section  III  indicates  areas  for  further  investigation.  The  approach 
followed  in  the  first  two  sections  is  to  generalize  or  extend  the  power 
spectral  density  function  for  stationary  stochastic  processes  so  that 
it  may  be  applied  to  non-stationary  stochastic  processes.  That  is,  the 
approach  follows  from  a  desire  to  obtain,  from  the  observation  of  a 
time-dependent  non-stationary  signal,  the  energy  contributions  at  one 
or  more  average  frequencies.  The  time  interval  of  observation  is  to 
be  sufficiently  large  that  an  energy  contribution  to  an  average  fre¬ 
quency  may  be  extracted,  but  not  so  large  that  the  process  changes 
appreciably  during  the  time  of  observation.  Thus,  through  successive 
observations,  the  change  in  average  frequency  may  be  observed. 

Consider  a  random  process  {K(t)}  consisting  of  an  ensemble  of  member 

functions  y,(t),  m  =  1 , 2 . M,  defined  for  the  interval  0  <  t  <  27’. 

For  this  process  the  autocovariance  '/'(t)  is  redefined  as 

DEFINITION  1 

where 


T 

+  ~  j  yjt  +  .  0  <  T  <  T 

0 


T 


y(t)dt 


1 


and 


r 

^  I  yjt)yjt  +  T)dt  ,  0<r<T  . 

0 

The  power  spectral  density  function  ^((o)  of  the  process  {^(t)}  is 
given  with  respect  to  '/'(t)  in  the  usual  way  except  that  the  interval 
over  which  the  Fourier  cosine  transform'  is  performed  is  finite,  hence 

DEFINITION  2 


^(ncoj.) 


cos  nii)j.TdT 


where 


<yjj.T  '  2n 

The  preceding  definitions  are  applied  to  a  subclass  of  non-stationary 
stochastic  processes  represented  by  those  processes  {F(t)}  consisting  of 
member  functions  "t  •=  1,2 . M  and  te[0,  2T]  which  may  be  rep¬ 

resented  by 


y,(t)  =  +  /(t) 

where  6^(t)  is  a  member  function  of  a  stationary  stochastic  process 
{0(t)}  and  /(t)  is  a  deterministic  function  inducing  non-stationarity 
and  is  statistically  independent  of  In  particular,  the  ensemble 

{F(t)}  is  said  to  satisfy 

Condition  lA] 

(1)  has  c.d.f.^  G(&^}  with  -  ro<  0^  <  -fee  for  every  m  »  1,2,  .  .  ,  ,  M 

and  each  te[0,  2T]  . 


f 


Cunqlatite  distribution  function  is  abbreTisted  c»defs 


2 


,(2)  E{0Jt)]  =  +  r)]  »  E[0]  »  J  ddG(0)  <  “  for  every  m  =■ 

1,2,  M  and  each  te[0,  27]. 

(3)  f{t)eL^  where  Lj  is  the  clas.s  of  square  integrable  functions. 

The  following  theorems  have  been  proven  in  Technical  Report  4. 
Theorem  1 

Given  an  ensemble  of  member  functions  y„(t)  “  0^(t)  +  /(f)  satis¬ 
fying  Condition  [4],  then  the  autocovariance,  >/'(t),  of  the  process  is 
separable  into  two  components  such  that 

'P(-r)  "  'Pg('r)  +  'pf(T) 


where 


^g{r)  -  E[0{t)0(t  +  T)]  -  £H0(f)] 

and 


'Pf('r)  =  -4^,  “  ^/(  «)'^/(  (+T) 

Theorem  1  shows  that  if  f(t)  -  0,  the  autocovariance,  '/'(t)  =  'pg(T), 
is  the  accepted  definition  for  stationary  stochastic  process.  Similarly, 
if  0^(t)=  0  for  all  m  and  if  /(f)  is  periodic  with  fundamental  period  T, 
then  the  autocovariance,  'P(t)  »  i/>^(t),  is  the  accepted  definition  for 
periodic  functions, 

Theorem  2 

The  power  spectral  density  function,  ,  for  an  ensemble  of 

member  functions  y„(f)  satisfying  condition  [d]  is  separable  into  two 
components  such  that 

I 

4>(nco^)  •  <^g(rvi>j.)  +  $^(nC0j,) 


3 


where 


and 


<t>a(nw^) 


cos  ncOjXdT 


T 

j  'ljf(T)  cos  nai^Tdr 
0 


This  theorem  follows  directly  from  Theorem  1  and  from  the  fact  that  a 
Fourier  transform  is  a  linear  operation. 

If  in  the  preceding  theorem  =  0  —  e.g.,  fit)  =  0 — then 


lim  0(ri<i>.)  »  lim  4>fl(nO)  )  «  ^oico) 

'  7*-40d  w  I  u 


TT 


^gir) 


cos  ctxr dT 


0 


ia  an  acceptable  definition  for  the  power  spectral  density  function  of 
a  stationary  stochastic  process;  and  if  'Pgi'^)  -  0  and  fit)  is  periodic 
with  fundamental  period  T,  then  ^(ncOj.)  «  <t?^inco^)  represents  a  definition 
for  the  power  spectral  density  function  of  a  periodic  function. 
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II  EXTENSIONS 


Generally,  in  applications  of  power  spectral  density  methods,  an 
additional  condition  is  imposed  on  ensembles  {l'(t)}  satisfying 
Condition  [,4]  so  that  estimates  of  $(n<Uj,)  may  be  obtained  from  a  single 
member  function  The  additional  condition  is  one  of  ergodicity 

and  is  stated  as  follows; 


Condition  [B] 

For  the  random  component  0^(t)  of  Condition  [A] 


E[e]  -  I  ddgiG)  =  .  lim|i  s  =  limir 


r-m  \M  «=1 


T 


ejt)dt 


E[0]  »  Ele  U+r)]  =  lim  i  2  0  (t  +  r)  »  limJ-T  ^  (t  +  rWt  I 

•  w  .=1  ■  r-®  It  J  ■■  r 

CO 

E[0{t)0(t  +  t)]  •  I  0{t)0(t  +  T)dHl0{t) ,  0{t  r)]  =  E[0^{t)0^{t  +  t)1 


limJ—  Z  0  (t)0  (t  +t)>  =  limJ~ 
"  -  T  7-coV 


0{t)0(t  +  T)dt 


where  the  time  averages  hold  for  each  value  of  m  =  1,2,  M,  and 

H[0(t),  0(t  +  t)]  is  the  joint  c.d.f.  of  0(t)  and  0{.t  +  t) 

Applications  assuming  ergodicity  imply  that  the  frequency  analysis 
is  derived  from  a  single  member  function,  indicating  that  an  estimate 
for  t//(r)  must  be  used  in  place  of  Definition  1.  Thus  the  following 
estimates  are  defined: 
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DEFINITION  3 


r 

B[0,(t)]  =  (  0,{t)dt 

^  -^0 

1 

(t  +  r)]  =  I  e  it  +  T)dt 

’  T  Jg  " 

.  T 

If 

E{6  (t)6  (t  +  t)]  =  ~  0  (^)0  (t  +  r)dt 

■  a  ^  J  »  ■ 

0 

A  A  A 

and  E[6^(t  +  t)]  are  estimates  for  EW],  and  E[6^(t)6^(t  +  r)} 

is  an  estimate  for  E[9  (t)6  (t  +  r)]  in  the  sense  of  Condition  [B] .  The  ex- 

*  .  •  .  A  . 

pressions  in  Definition  3  permit  an  estimate  for  (/’(t)  —  namely,  i/<(t) — which 

is  given  as  follows: 

DEFINITION  4 

A 

^  r  «+.T  )  »+t) 

and  for  the  random  component 
A 

A  ,  A  , 

Note  that  At)  =  ,{t)  .  Using  the  expression  for  ip  (t),  an  estimate  for 

’  ’  A  .  .  .  " .  . 

the  power-spectral-density,  <t(n&ij.) ,  is  given  by  Definition  5. 

DEFINITION  5 


A 


1. 

T 


(t)  cos  ncOjTdr 


Because  of  Conditions  [A]  and  [B]  and  Definitions  1  through  5,  $^(nct)j,)  is 
an  estimate  for  4>(nctij,)  in  the  sense  of  the  following  theorem. 
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Theorem  3 


Let 

{TCt)}  containing 


*  ?  Ss ,  , 

■  «. 5. 

member  functions  y 


Then  for  a  non-stationary  process 
^(t)  satisfying  Conditions  [A]  and  [B]., 


lim  {$„(nct)_)} 


[  riimj—  I 

I  [m-®  lu 


cos  ncOjXdT 


and  by  an  additional  limiting  process 


lim  lim  =  ®,(na)_)  +  4>fl(na)-) 

y*»CD  T-^CO  mi  j  I  CE  .  / 


Proof; 

By  Definitions  4  and  5 


T 

L  '0 


y.(*)y.(t  +  r)dt  “  — 


y^it)dt  y^(t+T)dt 


cos  nco^dr 


T  r'  T 


-I  ejt)e^t  +T)dt  +  -j  /(t)/(t  +T)dt  +- 
0  0 


+  r)clt 


&(t  +r)f(t)dt  --  f  &(t)dt 

t2  j  ■ 


&  (t  ■'  T)dt 


T2 

T2 


f(t)dt 

f(t  +  T)dt  — 

) 

7 

T 

f 

T 

6J^t  +  T)dt 

f{t)dt 

i 

) 

e  {t)dt 


o^dr 


f{t  +  T)dt 


7 


T 


“ l//(  «+t) 


cos  noifdr 


f\ 

I'Pg  (t)  +  flfy,)]  COS  mifdr 


where 


T  j  r 

^(y.)  '  r|  +  T)dt  — ^  e^(t)dt 


f(t  +  T)dt 


1 

+  -r 

T 


1 

-;7 


d,{t  +  T)f{t)dt - 1  6>.(t  +  T)rft 

?’■ 


Now  R(y  )  has  the  following  property  as  a  result  of  Condition  [B]  :■ 


fl  * 

'1  • 

)-  X 

|W  .«i 

? 

L  0 

Riy  )  cos  nco^rdr 


( litn 

U-® 


1  ^ 

-  2  R{yJ 
M  m’l 


COS  nco^rdr 


r 

Eie) 

0 

T 

0 

T  T 

E(e) 

fit  +T)dt  — Y 


Ei0)  f 

/(t+T)dt+‘-^  fit)dt 
0 


cos  nco^rdr  =  0 

Summing  $^(naij,)  for  m  =  1,2,  ...,  M,  and  taking  the  limit  as  W  ®,  we 
establish  the  first  part  of  the  theorem — i.e.. 


E[e] 


fit)dt 


ii™  i^Anci^-)}  “  li'"'  —  ^ 

u-u)  *  '  If-*®  [W  «•  1  J 


T  ^ 

*/'""'>*  M  [“:(« J.  ^*.‘^>}] 

A 


cos  nctijXdr 


Now  to  consider  the  limit  as  T  we  note  that 


lift)  (%  (t-)}  =  limjr  [  +T)dt  -  —  [  (9(t)dt 

f-4«  a  T'*^  j  7^  j 


0^(t  + 


From  Definition  3  this  is  equivalent  to 

lim  4g  (T)}  .  lim  a[e  (t)0  (t  +  T)]  -  El0(t)]heAt  +  T)]} 


T-*« 


T-'OD' 


and  by  Condition  [B]  and  Theorem  1, 


lim  {'Pg  (t)}  =  E[eit)e{t  +  T)]  -  EH0]  =  <Pg{r)  . 

T*-*® 


t^g  (t)} 

■ 

= 

5  limit 

as  1 

fl  * 

1 

IM  .=  1 

T 

lim  {pg  (t)  }  cos  nco  rdr 
r-®  ■ 


~  lim  « 

1 

JV-*® 

M  .=  1 

T 

_  fl 

^gir)  cos  nco^rdr 


^gincOj.)  as  a  result  of  Theorem  2 


Thus  lim  lim  {4>„(na)_)}  =  0,(n(y_)  +  ^g{nco-)  and  the  proof  of  Theorem  3 

If-*®  r— ro  mi  II  <7  1 

is  concluded. 


9 


As  given  by  Theorem  2,  4>^(rt<Uj,)  ijras  defined  with  respect  to  the  time 
of  observation,  T;  and  the  calculated  frequencies  are  multiples 
(harmonics)  of  the  fundamental  frequency,  However,  suppose  {T(t)} 

is  stationary  where  /(t)  is  a  periodic  function  with  fundamental  period 
^0  "  It  is  then  desirable  to  know  the  relation  between  the  calcu¬ 

lated  frequencies  and  the  harmonics  of  O}^.  Since 


'I 


<ty(nco^)  '  I  ‘=°®  noi^Tdr 


then 


1  cos-^->,rdr 


will  exhibit  the  fundamental  frequency,  co^,'  when  n  •>  N,  and  will  exhibit 
harmonics  or  multiples  of  o),  when  n  «  2N,  ZN,  4N .  The  notation  \p 


is  shorthand  for 


fH 


■  .  1  r"' 


”  ~  /(*)/(*  +  T)dt  - 


wher« 


f  = 


Wo 


UK 


f(t)dt 


NT, 


/(t  +  T)dt 


The  following  theorem  establishes  a  relation  between  the  calculated  and 
actual  frequencies. 

Theorem  4 

For  /(t)  e  Lj,  aperiodic  function  with  fundamental  angular  frequency 
o>0  and  fundamental  period  Tj, 


lim  )} 

T-HT.  ^  ^ 


n 

f\N 


10 


and 


(5N  +  2/3) 

{N+0)‘ 


where 

T 

Proof; 

Take 


(AT  +  P)To,  N  >  1.  Q  <  /3  <  1,  and  -  max_  ']/{t)| 


0<  K  To 


I  /(  tffit  +  T)c/t 


7 

--f 

7.2  I 


/(t  +  r)dt 


We  may  write  each  integral  in  two  parts,  as  follows; 


HT. 


^f(x)  •  - 


+  T)dt  f  +  T)dt 

»  k 


1 

T* 


#n. 


Lt» 


f(t)dt  +  [  /(t)dt 

wr„ 


r-  wr, 


0  r  ■ 

/(t  +  T)<it  +  f  /(t  +  T)dt 


HT„ 


- ;;  -  /(*)/(*  +  T)dt - + - /*  + —  f  +  T)dt 

S  +  PNT.  }  N  +  /3  N+/S  T  j 


_1_ 

7-2 


/vr,/+  f(t)dt 


AT./  +  I  f{t  +  r)dt 


11 


+  T)dt 


N  ,  /SAT  ,  1 

.4, 


(N+0)^  T 


«r„ 


f 

L  HT^ 


[fit)  +  fit  +  r)]dt  + 


fit)dt 


fit  +  T)dt\ 


If  this  expression  is  substituted  for  t/iy(r)  in  the  definition  of  ^^(ncOj.) 
then 


Vr)cosna,/rfr 


fit)fit+r)dt 


r 

-^[/VTo/j  [/(t+T) +/(t)]df 

/tT„ 


fit)dt 


f(t  +T)dt|^cosna>y.Tc^r 


Thus,  since  lim  /?  *  0,  taking  the  limit  as  T  NTf.  establishes  the  first 
T'^  ItT  ^ 

part  of  the  theSrem: 


lim  {<I',(nco  )}  «  "o) 


Next,  consider  upper  bounds  for  each  of  the  three  last  terms  in  the 
expression  for  4>^(na)j.); 


N  f  ip.j^  (t)  cos  noi^rdr 


iN 


jV _ f  1 

+  /?)T  J 

ItT. 


N  c  Ur 

(N+y3)rj  kJ 


1  r  fit)fit+T)dt-f^ 


cos  ncojrdr 


1  r 

(N  +  /S)^  ^0  L 


{N+l3)Tg 

( 

""0 

cos  ^  r 

J  * 

/V  +  >3 

"0 

\ctr  < 


20N  . 


iN  +  /9)2 


n  ' 


12 


-fl— 


cos  rtii^dr 


W  +  0)% 


(/*+y8)r, 


cos 


n 

Af+/S 


oi^rdr 


and  finally 


/(t)/(t  +T) 


J. 

fl 


T 

T 

T 

."j 

Wi/ 

[/(t+T)+/{t)]dt  + 

mdt  j 

f(t  +T)dt 

vcos  m>fdr 

/ 

»rg  / 

VTg  i 

»Tj 

J 

(N+/3)T, 


w)Tor 


(/V  +  /3)r, 


{n*/3)T^ 

I  + 


(Af*;9)rj 


(/V+/3y'T„ 


2V* 


(iv  +  /3)*r2 


(/*+/9)rg 

(W+/3)rg  1 

f 

n 

f„dt 

/»* 

cos  CO^T 

/TTg 

fffg  J 

rfr 


/S  „  .  ml  p^fl 


f  + 

N  +  ^  "  (/y  +  ^)2  +  1^7. 


{3N/3+fi^) 
{N  +  /?)^ 


Thus 


/  n  \ 

d>^(na)_)  - - <!>,( - o). 

^  (AT  +  /3)^  '■*■  P  I 


<  f 


,  0{5N  +  2P) 
(N  +  /3)* 


and  Theorem  4  is  proved. 

A  corollary  to  Theorem  4  gives  the  relations  between  i//^(t)  and 

N  =  1,  and  between  d>^(na)j.)  and  ^^Cnoig)  where  0  <  7’  £  Tg.  For  this 
investigation  set  T  »  yTg  or  coj.  »  fUg/y  where  0  <  y  <  1.  Remembering  that 


13 


(7-) 


+  T)dt  - 


where 


f 


+  r)dt 


then  the  corollary  may  be  stated  as  follows,. 


Carol lary : 


For  /(t)  £  1^2'  ®  periodic  function  with  fundamental  frequency  and 
period  ,  then 


lira  =  <l)^(no)j) 

T-Tq 


and 


<  5/^(7  “  i)/r^ 


where 


-  max 

/CO 

-  max 

/(t  +  t) 

0<  t<  Tjj 

“5OI-0 

and  0  <  T  < 


Proof ; 


4>f(T) 


T 

+  r)dt  - 

0 


r 


I  f(t)dt 
0 


T 

I  /(t  +  T)c/f 
0 


/(t)/(t  +  T)*  + 


TTo 


+  T)dt 


1 

■  ^0  H 

[  /(t)c(tf[  /(t)c(t 

-  ’■o 

/(t  +r)dt  +  j 

■yrj 

/(t  +T)dt 

J  J 

Lo  J 

-0  ■; 

'•o  -J 
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■yT. 


1  ,  y  -  1 

Ijfl  (t)  + 

V  „,2 


V'n  ~  ".  ~  /(*)/(‘  +  7-)rft 


y-  ^ 


yj’n 


7’n/ 


[/(t  +T)  +  f(t)]dt  + 


f(t)dt 


f(t  +  T)dt 


Now  since 


=  — 


rr„ 


cos  ncjo^rdr 


r^'n 


l/l.iT)  cos - Tc/t 

/  y 


substituting  the  preceding  expression  for  <p )  into  (p^CncCj.)  gives 


yr. 


■yj'n. 


r-r'  yr^  J 


y  ~  1  ,2  .  1 


y 


+  -  ,  — -  p  + 


y 


yT'o  1 


f(t)f(t  +  T)dt 


yr. 


rTo  YT^ 


y  ’’o 


[/(t  +t)  +/(t)]dt  -- 


f(t)dt 


f(t  +r)dt  >  cos - 

J  y 


yT„ 


y 


1  /'» 

2 


cos -  rdr 

y 


7r„  rr. 


{yT  J  J 

\  '  0 '  0  T 


/(t)/(t  +  T)dt  cos  -  rdr 

y 


3r2 


yJ7': 


0  0  r„ 


[/(t  +  t)  +  f{t)]dt  cos  - rdv-  + 

y 


rdr 
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and 


since 
Next , 
for  $ 

1 

■ 


1 


1 


(yT,y  -0 


r-  yy.n 


L  T. 


yTn 


fit  )dt 


f  {;t  +  T  )dt 


y 


■  rdr 


lim  (bfinco  ) 
T-T^  ’ 


lira  y  =  1„  and  the  first  part  of  the  corollary  is  established, 

r-r.. 

consider  upper  bounds  for  the  individual  terms'  in  the  expression 


yr. 


cos -  rdr 

’  y 


—  f 

ytf  J  L/( 

^  0  Tn 


f(t)f(t  +T)dt  -f^ 


cos -  Tdr 

y 


yX 


^0 

noj^ 

cos  ■ - T 

y 

rr. 

dT<^^  2fl 
T 


yTg  rT„ 


nCOn 


f{t)f(t  +T)dt  cos - Tdr 

y 


7T. 


a-y)T,fl 


cos - r 

r 


dr 


r(i  -  r)  ,2 


f 

J  u  t 


0  r„ 


[/(t  +t)  +  f(t)]dt  cos - rdr 

y 


2/^(1  -  r) 


7^7’ 


0  0 


rTo 

nco„ 

0 

cos - r 

7 

dr 


(1  -7) 


2/ 


2 

M 
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and  for  the  final  term 


y^o 

r 

yTo  y^o  -1 

1  f 

f  f 

ncUfl 

fit)dt  fit  +  T)dt 

cos - rdr 

(yTc)'  { 

J 

'•o  ’’o  J 

y 

(1  -  r)V^ 


nOfl 

cos  -  T 

0 

y 

dr  < 


(1  -  y) 


y‘ 


2 


The  proof  of  the  corollary  is  completed  since 


■—  "o) 


2(1  -  y)  ^r(i  ~  y)  ,2  ^  i  -  y)  ,»  ,(i  -y)^ 


■fl*' 

y2  y2 


./2  +• 
■'  M 


fi 


y^ 


r' 


2  liJlzL 1  -  .  C.2  (1  -y) 


<fl— — ^(4+r  +  i-r)  »  5/^ 

y2  y2 


The  power-spectral-density  techniques  are  applicable  to  non¬ 
stationary  processes  consisting  of  member  functions  which  may  be  reduced 
to  the  form  y^(t)  =  ^,(t)  +  f(t).  For  example,  consider  a  non-stationary 
process  {X(t)}  consisting  of  member  functions 


*.(‘)  =  >^.(t)g(0 


where  satisfies  Condition  [B]  and  Parts  1  and  2  of  Condition  [4] 

and  g(t)  satisfies  Part  3  of  Condition  [4],  Then  by  setting 


and 


y^(t)  «  In  xjt) 

djt)  .  In  KJt) 

fit)  ‘  In  git) 
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the  frequency  analysis  of  the  product  x^(t)  =  X^(t)g(t)  may  be  treated 
as  the  analysis  of  the  sum 

In  X.Jt)  +  In  g(t)  =  In  x^(t)  =  y,(t)  “  +  f(t) 

providing  the  dominant  frequencies  of  In  x^(t)  are  the  same  as  those  for 
This  restriction  is  stated  mathematically  for  functions 
observed  for  0  <  t  <  2T,  satisfying  the  following  theorem. 

Theobem  5 

Given  *„(*)  expandable  in  a  Fourier  series  as 


X  (t)  »  -  +  E  a  cos  nco^t  +  6  sin  na)-t 

,  2  n=l  "  T  n  T 


2 


[1  +  z_(0] 


where  the  and  6^  are  random  coefficients,  Oj  >  0,  co^T  ■=  27T,  and 
|*^(t)|  <<  1,  then  the  dominant  frequencies  contained,  in  the  set  ncOj.  for 
n  =  1,2,  ...  are  in  the  same  order  for  *,(T)  and  In  x^(t). 


Proof; 


Set  In  x^(t)  “  In  “  +  In  [l  + 


where 


z  (t)  = 


E  a  cos  noo  t  +  b  sin  nw-t 


n  ®  1 


T  n 


Since  lz^(t)|  <  1  and  >  0,  then  In  [l  +  z^(t)]  may  be  expanded  in 
powers  of  z^(t)  and  therefore 


“o  00  (-1)""^ 

In  X  (t)  =  In - +  E  -  z"{t) 

•  2  n=l  n  ■ 


The  condition  |z^(t)|  «  1  is  used  to  imply  that  the  first  term  of  the 
power  series  is  the  dominant  term: 
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.  0 

in  *,(t)  ”  in-^+z.(t) 

.  *^0  2  .  ® 

=  In  —  + - [  2  o  cos  nco-t  +  b  sin  nco_t J 

2  Og  „-i  "  ^  ^ 

which  provides  the  same  ordering  for  the  magnitudes  of  the  +  6*,  n  ^  1, 
for  in  *^(t)  as  for 

Thus  the  analysis  of  x ^{t)  =  X.^(t)g(t)  may  be  performed  on 
y^(t)  »  in  *^(t)  ”  ^  /(*)  by  the  usual  methods  of  taking  the 

Fourier  cosine  transform  of  the  autocovariance  of  y,(t)-  This  is  true 
because  the  choice  of  =  2tt/T  for  the  fundamental  frequency  of  g(t) 
carries  over  into  the  analysis,  since  g(t)  and  /(t)  ■  In  g(t)  both  have 
the  same  fundamental  frequency. 

Note  that  the  condition  «  1  is  equivalent  to 


—  »  2 
2  11“  1 


[al  + 


That  is,  the  dc  component  of  the  signal  must  be  greater  than  the  maximum 
amplitude  of  the  sum  of  the  non-dc  component  of  the  signal  so  that  the 
logarithmic  transformation  produces  a  real  function  which  is  physically 
interpretable . 
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Ill  ADDITIONAL  THEOREMS 


Consider  a  periodic  function  /(t)  with  fundamental  frequency 
and  period  which  is  expandable  in  a  Fourier  series — i.e., 


/(t)  =  +  2.  a  cos  n<D.  t  +6  sin 

2  n=l  "  On  0 


The  frequency  spectrum  of  f(t)  may  be  obtained  from  the  product  of  the 
finite  Fourier  transform  with  the  finite  conjugate  Fourier  transform. 
That  is,  as  the  finite  Fourier  transform  is  defined  by 


T 


0 


and  as 


the  finite  conjugate  Fourier  transform  is  defined  by 


0 


the  frequency  spectrum  is  given  by 


|F(nc<Jg)h  =  F{no)^)F*{nco^) 


For  periodic  functions  /(tleZ,^,  |F(nCt)jj)|^  =  4>^(nC0(j)  when  is 

evaluated  for  fit)  observed  over  the  time  interval  0  £  f  Tq .  However, 
is  functionally  different  from  O(na)j),  and  other  methods  of 
calculating  the  frequency  spectrum  of  fit)  may  be  investigated  (differing 
from  the  power-spectral-density  function  described  in  the  previous 
sections).  It  is  hoped  that  these  methods  can  be  extended  to  non¬ 
stationary  stochastic  processes  {F(f)}  satisfying  Conditions  [d]  and 

[sJ. 
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Consider  a  theorem  which  parallels  Theorem  4.  That  is,  let  /(t) 
be  observed  for  0  <  t  ^  T  where  T  »  (JV  +  /3)Tq,.  N  ^  1  and  0  £  ^  1. 

Again,  define  fUj.T  =  2vr  and  let 


T 

FincOf)  “ 


and 


F*  ( no), 


r 

,)  -  1 1  /(tye-^^T* 


dt 


Further,  since 


F(nojg) 


■  rl 

‘o  1 


/(t)e‘"“6'dt 


then 


Similarly, 


/ifr. 


r./"  \  1  r  ,/  V  •  <"  *  1 

'IS"")  ■  s?;  J  ■"  ■ 


F*(~  a)„)  =  -  /(t)e  dt 

\JV  7  JVT  j  ■ 


and  the  following  theorem  is  stated  and  proved. 

Theorem  6 

Given  /(t)  periodic  with  fundamental  frequency  and  period  Tg 
and  given  f(t)eL^  for  0  ^  t  ^  T,  then 


-  1# 


/g(4iV  +  ;3) 
(iV  +  /3)* 
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where 


Max  {|/(t)|} 

o<»  <rj 


Proof: 

Substituting 

and 


e  =  cos  noi^t  +  i  sin  nco^t 


„  ~  I  n  (u  _  e 

c  T  -  cos  rico  t 


i  sin  nco  t 


into  F(nC0j.)F*(nCt}j,)  gives 
P(noi^)F*  (no)\j,) 


'1  1 

r  2 

r  .  T  ri 

r 

1' 

T  J 
( 

/(t)  cos  no)  tdtj  j 

1  fit)  sin  ncii^tdt 

0 

1 


+  /S)ro  J 


ncu.t 


/(O  cos^dt 


1 

)vr„ 


nw, 


AT  + 


t 

7-"J 


*  /sir. 


iVT„ 


na>of 


V«+I3)T 


1  r  T' 

(»  *  m, 

)'rr„ 
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r  I  f'”""'*  ”"••  .,1 

1.(TT7I7;}_  .  ■'"’"'■T77  1 


wr„ 


[(/V  t  «r7  I  '“’  "”  Fr7'] 


r  1  1- 


(iV 


TwJ___  '] 


/»?■« 


r  1  f"’’®  ""o‘  1 

iw  *  ml] 


1 


(N  +  yS)?-.  J 


/(t)  sin 


/»r„ 


'*‘^®  *-  rftl 
/V  +  /3  J 
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! 


By  Schwartz's  inequality 


(I)  = 


f’” 

(N.  +  /Sir,  J 


(Ar+/3)r. 


(yv 


(«f+^)ro 

twA 


W+/J)r 


(W 


-4_  f 

+  /S)To  J  L  iV  +  /Sj 

#r„ 


and 


(H)' 


(iV 


ttIK) 


iV  +  p 


dt 


(N 


-^r 

*  L. 


«+/J)r„ 


[/(t)] *dt 


(AT 


1  f  7  1 


Thus , 


(/)*  +(II)^< 


(N  +  /3)r,  I 

ATT, 


(Af+>3)T5 


f^U)dt 


(N 


J  w^^)ro 

+  P)T,  J 


ffT, 


<  r 

—  J  M 


N  +  P^ 


Since 


ncOn  t  I 

[cos  — ^  111. 

'  N  +  p' 


""0  *  I  <r  1 
sin  — -  "  ■^-  151, 


N  +  P 


and 


l/(f)l  £  /»  , 
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(ben 


r_j _ f"\,„  i.F  1  r'”''"  ill  1 

w  ♦  /8)r,  J,  ^  AT  t  |S  ‘'‘J  [(Artflr,  L  “• 


+  2 


(iV 


1 

T^rJ  '“’•‘"5T7''' 


(jy 


-^f 

+  /3)r,  J 

®  Afr, 


W+AjTj 


nCt)*  t 

/(t)  .in^rft 


<4/2 


/SiV 


(iV  +  = 


and 


^  yfl(4iV  +  /3) 

~  +  /3)®  *'* 


As  in  the  case  of  Theorem  4|  a  corollary  to  Theorem  6  exists.  The 
statement  of  the  corollary  follows  without  proof  since  the  proof  parallels 
that  for  the  corollary  to  Theorem  4. 


Corollary  1 

Given  /(t)  periodic  with  fundamental  frequency  and  period  T, 
and  given  f{t)eLj  for  0  <  t  <T,  then 


F(na>j)F*(nOij)  ”  ~ "o)  "c)  | 


A  (S  -  7)(l  -  y) 


where 

f,  -  Max^  {|/(e)l)  ,  T  -y^  , 


2S 


and 


0  <  y  <  1 

For  stochastic  processes,  consider  the  process  {l'(t)}  containing 
member  functions  y,(t)  •  ^„(t)  +  /(t)  satisfying  conditions  [A]  and 
[B]  as  before;  and  define  a  sample  function  g,(t)  as  follows. 

DEFINITION  6 

Let  a  sample  function  g^(t)  be  constructed  for  each  member  function 

y.(*)  by 

g.(0  -  +  /(t)  -  A,,,,] 

-  +  [/(t)  -  A^,,,] 

where 

T 

'‘/.(O  ”  T  1 

*^0 

1  ^ 

(t)  ’  r  I 

tr 

<1 

'  T j  • 

'^0 

'fh^  Fourier  transform  of  g,(f)  for  0  ^  t  ^  T  and  for  d^{t)€L2  and 
is  given  by 


26 


i 


T 

r,  y  V  If  .  .  in<0-T 


-ff  («.(<)  -'1«  ,.,1‘'"“''‘'<  *\\  UM 


where  oijT  =  27f;  and,  similarly,  the  conjugate  transform  is  given  by 


/  \  ""t'  , 

g^( t)e  ^  dt 


■  ■?[  w.ci  I-"""'’*  ■ 

0  *  "^0 

The  expected  value  of  the  power  spectrum  is  given  by  the  mathematical 
expectation  of  the  ensemble  of  the  products  of  the  Fourier  transform  of 
each  stochastic  member  function  with  its  conjugate  transform— i . e . , 

£[lf(n)|*]  =  £[F.(n)f;(n)] 

2  ,  T 


«  « 

T 


®[{t  j  nco^trft}  -  I  cos  ncOj.td^- 

0  0 

T  T  2_ 

'  If  I  n<o^td^  {7  j  nco^td^  1 

0  0 

r  r  2  T 

+  ^  J  nct>j.td(j  *  2 f  (f) 

'^0  0 

1  r  1  ^  h 

•  {7;  J  =i"  noy^tdt^  +  j  sin  I 

0  0  * 

T  2  T  2 

"A^(,,]cosna)j.tdtj^  +  1/(0  sin  nOYt*} 
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Since  £[5^(t)]  »  ElAg  (,p.  the  precedi  ng  expression  reduces  to 


_  7*  2  7* 

£t|f(n)|*]  «  E  IJ— I  0^(t)  cos  nOJj.tdtj-  j  “  B  1“  J  cos  «aytdt| 


r  1  r  2^  r 

^  ^  11  "^T**^*}  J  ~  [  si"  nco^tdtj 

0  4 

r  2 

*  Ir  /  ""r**} 

0 

T  2 

^  {t  J  [/(f)  sin  ni)j.tdt| 


Next,  consider 


«  ^  [M  f^.(f  )0,(T)  [cos  rxu.t]  [cos  na>,T]dtrfr 

0  ^  0  0 


E[6^(t)]  '  •^— Jcos  nw. 

0 


and,  similarly, 


^  [(t-  j  ‘  ^  ®i"  ] 


Vtdf} 


B[0*(t)]  •  II  J  sin  noij.tdt^ 
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1 


Thus , 


r  r  2  T  2-1  p 

£[l/(n)|^]  “  [{7  J  j  nc^>j.td^  J  •  “  £^t^] 


(7]  "'^/(.t)^“®  "“r‘* 

0 


} 


T 

{■^  j  [/(«)  sinnwj,ta[t| 


Finally,  since  ci)^T  =  27t 


I 

T 


cos  noiftdt  = 


1  for  n  =  0 
0  for  n  >  1 


and 


sin  ncOftdt  =  0  for  all  n 


and  the  following  theorem  has  been  proved. 

Theorem  7 

For  the  sample  functions  of  Definition  6,  the  expected  value  of 
the  power  spectrum  is  separable  into  a  stochastic  component  and  a 
deterministic  component  —  i.e., 


|F(n)P  =  |.FJ*  +  |F^(n)|* 


where 


E[0^l  -  for  n  =  0 

0  for  n  >  1 
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and 


1  r^'  1 

“  [/{t)  ~  nco^tdtj 


2 


sin  n^jtdt 


The  result  stated  in  this  theorem  is  to  be  expected  since  /(t)  and  the 
&^(t)  are  statistically  independent. 

In  the  following  discussion,  possible  extensions  for  the  use  of 
transform  methods  are  considered  for  the  analysis  of  non-stationary 
deterministic  functions  fit)e.L^.  For  example,  let  /(t),  0  <,  t  £  Tj, 
be  defined  by 


where 


/(f) 


<w(f  ) 


cos  [ct)(t)t]  = 


Clu  cos 

\ 

*^0^0 

“  27T 

Define  a  transform  by 


and  a  conjugate  transform  by 


F*(n) 


Substituting 
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in  F(n)  and  F* (n) gives 


Fin)  .  — 


L^o 


and 


F*in) 


J  n 


0  ''0 


"  r  I  " 

‘  0  Jo 


Completing  the  square  in  each  exponent. 


To  T, 

—  ■it  =  -  +  t+ - 

To  r„  4  4 


+  !5.  V  - 

'  2  y)  “4 


then 


±<(n  +  l)a)j 


,  <  I^V  ^0 

2  j  ±Mn+l>a.(,— 


ii(n-l  )<i). 


f(n)  =»  — 


"t(n+l)^/2 


'V^(n  +  1)  J 


e  de  + 


-i(n-l)Tr/2 


y277-(n  -  1)  J 


e'^  dB 


and 


F*(n) 


,  i(n+l  )7r/2 


i2-n(n  +  1)  J 


e"‘®  d0  +■ 


y27r(n  -1)  J 


e"‘®^  dB 


dt 
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where 


and 

Thus, 

F(n)F*M 


Finally,  it 

FMF*(n)  = 


a  »  —  ✓27r(n  +  1) 

,  3  - 

6  =  —  /277(n  +  1) 

c  ■  y  /27r(n  -  1) 


d  -  —  y2iT(n  ~  1) 


27r(n  +  1) 

V. 


[cos  6^  +  i  sin 


[cos  -  i  sin  d^]dd 


+ 


2n/n^  -  1 


+ 


1 

27r(n  “  1) 


d 

J  [cos  +  i  sin  d'^'idO 


in  be  shown  that 


^  [  L/TTT 


e^de\ 


L/TH 


sin  e^de  -  ^  [  sin  ji-  ,  n  >  1 

/n  -  1 
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That  is,  this  extension  to  a  transform  permits  the  calculation  of 
|f(n)M  for  the  function  f(t)  -  coa  coJ.(t  +  in  terms  of 
Fresnel  integrals. f 

To  see  that  F(n)  may  be  interpreted  in  terms  of  a  Fourier  trans¬ 
form  consider 


where  CJ^  is  a  constant.  Then,  we  may  write 

r  f 

ico.  Iff  n  .  in^^t 

f(t)e  ^  e  ^  dt  =  y*j  /(f)[cos  a)^t^  +  i  sin  a>jt*]r  '  dt 

0  0 

By  setting 


*(  f  )  =  /(  t  )  cos  * 
y(t)  =  f(t)  sin  ojjt* 


and 

t(t)  »  *(()  +  iy(t)  , 

then 

T 

F(<a)  ^  z{t)e  ^  dt 

0 


and  F{.o})  is  the  Fourier  transform  of  the  complex  function^  z(t). 

The  preceding  analysis  suggests  considering  finite  transformations 
of  the  form 

=  -J  /(£)e  ®  dt 

0 


*For  n  -  I, 


e  -  H  and  two  of  the  integral*  in  F(!)F*(l)  are  lero* 


For  n  =  0,  F(0)f(0)  = 


/(t)rff) 


3.3 


where 


P  M 

P(t)  =  t  +  2  . 

P  =  2 

The  ct)^  are  constants  to  be  adjusted  to  find  a  "best  fit"  with  non¬ 
stationary  functions  fit),  0  <  t  ^  T,  for  a  series  similar  to  the  Fourier 
series.  As  a  first  step  in  this  direction,  the  function  Fin)  is  shown 
to  be  the  Fourier  transform  of  the  complex  function  r(t)  =  xit)  +  iyit) 
where 


P 

'  /  P  \  1 

{  P  \ 

zit)  •  fit)  exp 

i  1.  CJ  tP 

p  =  2 

-  fit) 

cos  1  Z  o)  t'’;  +  sin 
\(.=2  /  ' 

^P=2  "  / 

=  xit)  +  iy(t) 

For  0  ^  t  _<  T,'  the  Fourier  transform  is  represented  by 

FM  ■  ij  d, 

and  a  sufficient  condition  for  the  existence’  of  Fico)  is  that 

1 

-  |j(i)Ui  <  “  . 

Consider  the  modulus  of  z(f.), 


f 

/  P  \ 

/  P  \ 

U(t)l  = 

cos  (  Z  &)  f'’ 
\p=J  "  / 

1  +  i  sin  (  Z 

P 

Z  cu  t'’ 

p-2  P 


=  \fit)\  . 

Thus,  since  /(t)ef,j,  \zit)\eL^  and  the  transform  exists.  The  existence 
of  F*iM)  is  demonstrated  in  a  similar  fashion.’ 
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